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In order to effectively realize and control the critical coupling, a graphene-based hyperbolic metamaterial 
has been proposed to replace the absorbing thin film in the critically coupled resonance structure. Our 
calculations demonstrate that the critical coupling effect (near-perfect light absorption) can be achieved at 
the near- infrared wavelength by using this layered structure, while the critical coupling frequency can be 
tuned by varying the Fermi energy level of graphene sheets via electrostatic biasing. Moreover, we show that 
the critical coupling frequency can be tuned by changing the thickness of the dielectric or layer number of 
the graphene sheets in the unit cell of the graphene- dielectric HMM. The optimization performance has also 
been indicated, which may offer an opportunity towards the experimental designs of high efficient graphene 
based critical coupling devices. 



Graphene has attracted intense scientific interests owing to its outstanding optical properties including 
the strong light- graphene interaction, the broadband and high-speed operation, giant nonlinear Kerr 
index (several orders of magnituede larger than conventional bulk media) and the controllable 
optical properties due to tunability of the charge carrier density by applying the bias voltage upon gra- 
phene 1 " 8 . Graphene also exhibits wavelength insensitive absorption in the visible- and near-infrared frequen- 
cies, which makes graphene a potential candidate as a building block for optical detectors and photovoltaic 
devices. However, absorption of the single-layer graphene remains very weak. Some methods had been 
proposed to enhance its absorption. Thongrattanasiri et al. demonstrated that 100% light absorption could 
occur if a single patterned sheet of doped graphene is designed instead 9 , Ferreira et al. reported that the 
absorption of graphene can reach up to 100% within Fabry-Perot cavity sandwiched graphene structures 10 or 
in graphene-metamaterial corrugated structures 11 , Nikitin et al. showed that graphene periodic antidot arrays 
could provide a strong electromagnetic absorption at both microwave and THz regions 12 , Nefedov et al. 
demonstrated that 100% light absorption could also be achieved in a graphene-based asymmetric hyperbolic 
metamaterial (HMM) 13 . In this paper, we show that near-perfect light absorption can be obtained by the 
critical coupling within graphene-based HMMs. 

A critically coupled resonator (CCR) is a thin-film structure that can absorb nearly all of the incident 
electromagnetic radiation, leading to the null scattering. It was firstly reported by Tischler 1415 , who experi- 
mentally demonstrated that a 5 nm thick film of high oscillator strength J-aggregated dye separated from a 
dielectric Bragg reflector by a dielectric spacer can absorb more than 97% of the incident light at the 
wavelength X = 591 nm 14 . In the following, Gupta et al. replaced the polymer absorbing layer by a metal- 
dielectric composite film, and showed that the critical coupling (CC) frequency can be tuned by changing the 
volume fraction 16 . For a larger volume fraction, CC at two distinct frequencies can be obtained due to the 
strong coupling. Subsequently they extended their studies to oblique incidence for both the TE and TM 
polarizations 17 , Fabry-Perot cavity with metamaterial mirrors 18 , or layered media with Kerr nonlinearity 19 . 
However, a critically coupled system has limited flexibility for a given absorber material. The tunability of CC 
is intrinsically important for applications in optical devices. It is desired to tune CC in a fixed configuration 
or device by an external field. For this purpose, introducing graphene sheets into a CCR seems to be a good 
approach for achieving CC. Here, we replace the polymer absorbing layer or metal -dielectric composite film 
by the graphene-based HMMs 20 " 23 composed of stacked graphene sheets separated by thin dielectric layers. It 
is found that such a graphene-based HMM can support controllable CC at a single wavelength. 
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Figure 1 | (a) Geometry for a graphene-based hyperbolic metamaterial (HMM). (b) A critically coupled resonator with the graphene-based 
HMMs and perfect mirror spaced by the air. (c) The critically coupled resonator with the graphene-based HMM and the dielectric distributed 
Bragg reflector. 



Results 

Graphene-based hyperbolic metamaterials. Our graphene-based 
HMM is representatively shown in Fig. 1(a). It is a graphene 
multilayer structure in which separations between graphene sheets 
are assumed to be filled with the dielectric. For the HMM in air, as 
shown in Fig. 1(a), the absorption is smaller than 50% generally and 
the perfect absorber cannot be constructed. In order to realize the 
perfect absorption, we try to construct a CCR by introducing a 
lossless metallic mirror to the back of the HMMs, as shown in 
Fig. 1(b). In the experiment, lossless metallic mirrors may be 
replaced by more realistic dielectric Bragg reflectors (DBRs). The 
designed CCR has been shown in Fig. 1(c), which is consisted by a 
graphene-based HMM separated from a DBR by a dielectric spacer. 
The graphene-based HMM acts as a thin absorbing film, which is 
created by alternating layers of graphene with the surface 
conductivity a and slab thickness t g and conventional dielectric 
with a dielectric constant s d and slab thickness t d . The period of 
the unit cell of the layered structure is t = t d + t g . The graphene's 
surface conductivity can be calculated by adopting the Kubo 
formula 24 " 26 . Without considering the external magnetic field, the 
isotropic surface conductivity a of graphene can be written as the 
sum of the intra-band a intra and the inter-band term a in t e r> where 
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where co is the frequency of the incident light, E F is the Fermi energy, 
x is the electron-phonon relaxation time, and Tis a temperature in K. 
e, H and k B are the universal constants related to the electron charge, 
reduced Planck's constant, and Boltzmann constant, respectively. 

The Fermi energy E F = hv F (nn 2 D) 1 ^ 2 can be electrically controlled 
by an applied gate voltage due to the strong dependence of the carrier 
density n 2 r> on the gate voltage, where v F = 10 6 m/s is the Fermi 
velocity of electrons. Obviously, a is highly dependent on the work 
frequency and Fermi energy, which could provide an effective route 
to achieve an electrically controlled CC phenomenon. 

Under the assumption that the electronic band structure of a 
graphene sheet cannot be affected by the neighboring graphene 
sheets, graphene's effective permittivity s g can be written as 



= !+■ 



SQCOtg 



(3) 
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where s 0 is the permittivity in the vacuum. Considering the Cartesian 
coordinate systems shown in Fig. 1(a), layered structure are illumi- 
nated from the left with space- time dependence of fields as Aexp(ik z z 
+ ikxX — icot). In the sub- wavelength limit, graphene- dielectric 
layered structure can be treated as a homogeneous effective medium 
with anisotropic permittivity sensor, 



(4) 



where s x = (f g /s g +f d /s d ) \ s y = s z =f g *e g + fd*£d>f g = t g /tandf d = 
t d /t are the filling ratio of graphene sheet and dielectric, respectively. 
id 

If t g <C tj, s y = s z ~s d H , which is consistent with the results by 

averaging the effective displacement current (including both the dis- 
placement current in the dielectric slab and conduction current in the 
graphene sheet) over the associated electric field in a unit cell 20 ' 21 . 

For the TM waves propagating in the graphene-dielectric layered 
structure, the spatial dispersive curve is given by 




(5) 



where ko is the free- space wave vector, k x and k z are the wave vector 
along x and z direction in the graphene-dielectric layered structure, 
respectively. If e x e z < 0, the dispersive curve of Eq. (5) is hyperbolic 
type, and we call this graphene-dielectric layered structure as HMMs. 
Since the effective permittivity s x and s z depend on the Fermi energy, 
the dispersion characteristic of HMM could be tuned by changing the 
Fermi energy of the graphene sheets. 

As an example, in Fig. 2 (a) -(d) we plot the effective permittivity s x 
and s z as a function of the operation wavelength at the different 
Fermi energies at the near-infrared (NIR) frequency. Here, we 
assume that the Fermi-energy of graphene E F = 0.46, 0.48, and 
0.50 eV, respectively. The other parameters are T = 300 K, x = 
0.5 ps, t d = 8 nm, s d = 11.7; t g = M*0.35 nm, where M is the layer 
number of the graphene sheets, and M = 1 is for the single-layer 
graphene. It is seen from Fig. 2(a) and (b) that the real part of s z is 
positive, and both the real part and the imaginary part of s z vary with 



wavelength slowly. As the Fermi energy is increased, both the real 
part and imaginary part of s z decrease slightly. However, s x exhibits 
resonant behavior in NIR frequency. Re(s x ) < 0 and Im(s x ) become 
very large near the resonant wavelength ^ res . For E F = 0.50 eV, X res 
= 1.534 um, at this frequency, Re(s x ) < 0 and Re(s z ) > 0. Hence Eq. 
(5) denotes a dispersion curve of hyperboloid. Furthermore, we find 
that the resonant wavelength of the hyperbolic dispersion curve can 
be electrically controlled by an applied gate voltage on the graphene 
sheets. It is clear that the resonant behavior of s x can be tuned by 
varying the Fermi energy. Increasing Fermi energy, the resonant 
wavelength ^ res shifts to the smaller wavelength; while decreasing 
Fermi energy results in the red shift of the resonant wavelength X res . 

Critical coupling phenomena with graphene-based HMMs. The 

resonant behavior of s x can be used to design a CCR. We replace the 
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Figure 3 | The absorption of critically coupled resonance as a function of 
wavelength at different Fermi energy E F , where t s = 600 nm, t H MM = 
83.5 nm, and the other parameters have the same values as those in 
Fig. 2(c). The dotted lines are the corresponding absorption of the HMMs 
on the glass substrate. 
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Figure 4 | The magnetic field distributions in the critically coupled resonance at different Fermi energy E F , (b) is the enlarged view of (a) in the 
graphene-based HMMs and spacer, (c) is the magnetic field distributions gotten by FDTD method for E F = 0.50 eV and X = 1.534 um. E F = 0.50 eV 
and X = 1.534 um, and E F = 0.48 eV and X = 1.594 um are for the solid-line and dashed-dotted line, respectively. Other parameters are as in Fig. 2, the 
red dashed-dotted lines in (a) and (b) and solid-line in (c) show the interface of different materials. 



polymer absorbing layer or the metal- dielectric composite film in the 
conventional CCRs by the graphene-based HMM, as shown in 
Fig. 1(a). The DBR (AB) N A is constituted by two dielectrics A and 
B with permittivity s A and s B and thickness t A and t B respectively. In 
the numerical calculations, we choose the following parameters: N = 
10, s A = 5.7121 and s B = 2.6244. The thickness of dielectrics A and B 
are chosen to correspond to X/4 plate: t A = ^ c /(4n A ) and t A = XJ 
(4n B ), where the central wavelength is chosen as X c = 1550 nm, 
tia = V^A' an d Kb = y/ss. To easily design and adjust the spatial 
distance t s , we assume that the space between the graphene-based 
HMM and DBR is filled by air. 

Next, we show that CC can be realized in our structure. In Fig. 3 
we give the absorption of CCR as a function of wavelength at the 
different Fermi energy E F . The thickness of spacer layer t s = 600 nm 
and the thickness of HMM t H MM = 10* (tg + t d ) = 83.5 nm are 
chosen to achieve nearly 100% absorption at the central wavelength 
^ c = 1.534 um and E F = 0.5 eV, and the other parameters have the 



same values as those in Fig. 2. Compared to the absorption A ~ 40% 
of the HMMs on the glass substrate, the absorption from the CCR is 
enhanced about 2.5 times. Hence, our structure becomes very effi- 
cient to maximize light absorption. Another advantage of our struc- 
ture is the tunable CC frequency by changing the Fermi energy levels, 
as shown in Fig. 3, which are in consistent with the dependence of 
the effective permittivity s x on the Fermi energy as shown in 
Fig. 2(c). When we decrease the Fermi energy E F , the absorption 
peak moves to the longer wavelength side; when we increase the 
Fermi energy E F , the absorption peak moves to the longer wave- 
length side. For E F = 0.46 eV, the central wavelength is 1.661 (im; 
however for E F = 0.48 eV, the central wavelength X c is 1.594 |im. 
This property suggests that the frequency of absorption peak can be 
engineered by the Fermi energy of the graphene sheets. Although the 
absorption remains still high (>90%), either the decrease or increase 
in the Fermi energy level leads to the departure from the CC 
phenomenon. 
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Figure 5 | The contours of the absorption of critically coupled resonance as a function of thickness of spacer layer t s and thickness of HMMs 
at the different Fermi energy E F . Where E F = 0.52 eV and X = 1.478 um in (a), E F = 0.50 eV and X = 1.534 um in (b), and E F = 0.48 eV and 
X = 1.594 um in (c). Other parameters are as in Fig. 2. 



To understand the almost 100% absorption of the incident energy 
by the CCR at the CC frequency, we have plotted the magnetic field 
distributions inside the layered structure at the CC wavelength X = 
1.534 um for E F = 0.50 eV based on the transfer matrix method 
(TMM), as shown in the solid-line in Fig. 4(a) and (b). To confirm 
the correctness of TMM, in Fig. 4(c) we also provide the simulated 
magnetic field distributions by using the Finite-difference time- 
domain (FDTD)(The time evolution of magnetic field has been given 
in the Supplement materials). The magnetic field amplitude has been 
normalized to the incident magnetic field amplitude. One can see 
that the magnetic field is evanescent light wave in DBR due to the 
total reflection, however the magnetic field is strongly enhanced and 
localized in the graphene-based HMM and spacer. The large intens- 
ity enhancement and large imaginary parts of s x in the absorption 
layer of graphene-based HMM, explains the CC phenomenon (near- 
perfect light absorption) at the CC wavelength X c = 1.534 um for E F 
= 0.50 eV. Additionally, for comparison, in Fig. 4(a) and (b) we have 
plotted the normalized magnetic field distributions for X = 1.594 um 
at E F = 0.48 eV. It is clear that those magnetic fields are weaker than 
the field at E F = 0.50 eV, explaining why there is a departure from 
the CC phenomenon at E F = 0.48 eV. Hence, in order to achieve the 
CC at the different Fermi energy levels, we must optimize the phys- 
ical parameters, such as, the fill factors of the dielectric and graphene 
sheet in the HMMs, the thickness of spacer layer t s , and the thickness 
of HMMs t H MM- However, to realize the tunability of the CC fre- 
quency by changing the Fermi energy in a fixed configuration, we 
intend to finely adjust the distance between the graphene-based 
HMM and DBR on glass substrate, i.e. the thickness of spacer layer 
t s . The optimization is given as following. 

Recover the CC at the different Fermi energy levels. To recover the 
CC at the different Fermi energy, we have plotted the contours of the 
absorption of the CCR as a function of the thickness of the spacer 
layer t s and the thickness of the HMMs at the different Fermi energy 
E F , as shown in Fig. 5. The red region gives the absorption larger than 



99%, here we consider it as the CC region. It is obvious that both t s 
and t H MM should be decreased with the increase of Fermi energy. To 
easily adjust this property in the experiments, we intend to find the 
common area of the thickness t H MM for the CC at the different Fermi 
energies. It is found that this region is located near t H MM = 88 nm. 
The optimized results have been demonstrated in Fig. 6. Here we 
keep the thickness of HMM t H MM — 88 nm fixed, and adjust the 
distance t s . According to the optimized results in Fig. 5, we choose t s 
= 590 nm for E F = 0.50 eV, t s = 610 nm for E F = 0.48 eV, and t s = 
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Figure 6 | The absorption of critically coupled resonance as a function of 
work wavelength at different Fermi energy E F . Where t s = 590 nm for E F 

= 0.50 eV, t s = 610 nm for E F = 0.48 eV, and t s = 650 nm for E F = 
0.46 eV, t H MM is kept as constant 88 nm. Other parameters are as in Fig. 2. 
The gray dashed-dotted line gives the reflection of DBR on the glass 
substrate. 
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Figure 7 | The critical coupling phenomena at two distinct wavelengths at 
different Fermi energy E F . ForE F = 0.50 eVweusedt H MM = 16* (t g + t d ); 
for Ep = 0.48 eV we used t HM M = 17*(t g + t d ); for E F = 0.46 eV we used 
t H MM = 19* (t g + ta), where M = 1 and t s is kept as constant 500 nm. Other 
parameters are as in Fig. 3. 



650 nm for E F = 0.46 eV. It is clear from Fig. 2 that, for X = 
1.661 um at E F = 0.46 eV departure from CC, one can recover CC 
just by adjusting the distance t s = 650 nm. For X = 1.594 um at E F = 
0.48 eV, CC resonance can be recovered by increasing the distance t s 
to 610 nm. Therefore, our CCR has the potential to achieve tunable 
CC frequency in a fixed configuration. 

Critical coupling phenomena at two frequencies. From Fig. 3, it is 
found that there are two absorption peaks in the layered structure. 
In Fig. 6 we have recovered the CC at the different Fermi energy for 
one peak. But how can we realize the CC phenomena at two distinct 
frequencies? By optimizing the parameters of the structure, 
including t s and t H MM> we nave plotted the absorption coefficient 
as a function of wavelength, as shown in Fig. 7. We have used 
different sets of values for different Fermi energies. For E F = 
0.50 eV we used t HM M = 16* (tg + td); for e f = 0-48 eV we used 
t H MM = 17*(t g + td); for E F = 0.46 eV we used t HM M = 19*(t g + 
ta), where M = 1 and t s = 500 nm. It is clear that for optimal values 
of the thickness of graphene-based HMMs we can obtain CC at two 
distinct frequencies. 

Discussion 

In addition to Fermi energy, the resonant behavior of s x also depends 
on the fill factions of dielectric and graphene sheet. Decrease in 
thickness of dielectric t d or increase in layer number of graphene 
sheets M leads to the decrease of fill factor of dielectric f d and the 
increase of fill factor of dielectric f g , so that the graphene sheets are 
becoming more and more important in the HMMs. These properties 
are significant for controlling CC due to the strong dependence of the 
CC behavior on the resonant behavior of s x , while the results are 
shown in Fig. 8. The solid lines are the CC phenomenon at the central 
wavelength X = 1.534 um and E F = 0.5 eV, the parameters are the 
same as those in Fig. 3. In Fig. 8(a), it is found that the central 
frequency of the absorption peak shifts towards longer wavelength 
side with increasing the layer number of the graphene sheets M. The 
absorption peak departs from the CC phenomenon, which may be 
optimized by varying t H MM or t s as shown in Fig. SI of the 
Supplementary materials. In Fig. 8(b), we find that the central fre- 
quency of the absorption peak shifts to longer wavelength side with 
decreasing t d and shifts to shorter wavelength side with increasing t d , 
the optimization result is given in Fig. S2 of the Supplementary 
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Figure 8 | The absorption as a function of work wavelength at the 
different layer number of graphene sheets (a), thickness of dielectric in 
graphene-based HMM (b) and relaxing time (c). Where td = 8 nm and x 
= 500 fs in (a), M = 1 and x = 500 fs in (b), and t d = 8 nm and M = 1 in 
(c), other parameters are as in Fig. 3. 

materials. In Fig. 8(c), we show the absorption as a function of 
wavelength at different relaxing time x. It is clear that the central 
frequency of the absorption peak keeps fixed and the absorption 
decreases with the decreased or increased relaxing time, the optim- 
ization result is given in Fig. S3 of the Supplementary materials. 
According to above analyses, it is demonstrated that the dependence 
of absorption on the resonant behavior of s x provides more degree of 
freedom to control the CC phenomenon at the NIR frequencies by 
using graphene-based HMMs. 

In conclusion, we have presented a graphene-based HMM as the 
absorbing thin film to realize the CC phenomenon in the NIR fre- 
quency regime and discussed the controllable properties of this phe- 
nomenon in the CCR. It is found that the CC phenomenon depends 
on the resonant behavior of the effective permittivity of HMMs, and 
the CC frequency can be effectively tuned by changing the Fermi 
energy level applied on the graphene sheets. However, near-perfect 
absorption departs from the CCR if the Fermi energy is significantly 
changed. Hence, to recover the CC phenomenon optimization on the 
thickness of dielectric, the layer number of graphene sheets in 
HMMs, and the relaxation time of graphene sheet need to satisfy 
certain requirements. The optimization results are also presented. 
We believe that the controllable CC phenomenon and near-perfect 
absorption at the NIR frequencies (telecommunication band) could 
find potential applications in optical communications, optical detec- 
tors and photovoltaic. 

Methods 

The numerical calculations for the effective permittivities of the graphene-based 
hyperbolic metamaterials, the theoretical absorption spectrum, and the optimization 
results of critical coupling phenomena in the critically coupled resonant containing 
graphene-based hyperbolic metamaterials were generated using Matlab code based 
on the transfer matrix method. Let a plane wave be injected from vacuum into the 
critically coupled resonant at normal incidence with +z direction. Generally, the 
electric and magnetic fields at any two positions z and z + Az in the same layer can be 
related via a transfer matrix, 
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/ cos (k iz Az) —jqiz sin (k iz Az) \ 
l t (Az,co) I _j_L s i n (j Cig Az) cos (k iz Az) j 



where i denote the different media, k iz = J k^Si — k\ is the components of the wave 
vector along the z axis in the isotropic medium. For the TM polarization, 
k z = kl& x — (& x / & z )k 2 x and q z = k z /k 0 8 x in the hyperbolic metamaterial, q iz = kjko&i 
in the isotropic medium. The reflection and transmission coefficients can be calcu- 
lated as 

_ [X M +X h2 * pi] * pi — [X 2 ,i +X 2;2 * pi] 
~ [X M +Xi >2 * pi] * p\ + [X 2 ,i +X 2 , 2 * p L ] 



[Xi A +X lj2 * p L ] * p l + [X 2jl +X 2 , 2 * pi] 

Hence, the absorption can be obtained as, 

A=l-\r\ 2 -(p L / Pl )\t\ 2 , 

where Xy(/,j = 1,2) are the matrix elements of X ir j(co) = YlMi(Az,co), which 
represent the total transfer matrix connecting the fields at the incident end and the 

exit end; pi = \jk\z x -k\j ei ,p L = \Jk\z L -k\l &i, Si and s L are the permittivity of 
the media before the incident end and the space after the exit end, respectively. 
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